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Abstract 

There are many natural, physical, and biological systems that exhibit multiple time scales. 
For example, the dynamics of a population of ticks can be described in continuous time during 
their individual life cycle yet discrete time is used to describe the generation of offspring. These 
characteristics cause the population levels to be reset periodically. A similar phenomenon can 
be observed in a sociological college drinking model in which the population is reset by the 
incoming class each year, as described in the 2006 work of Camacho et al. With the latter as 
*rlopezl4(§asu.edu, ^ sks@asu.edu, ^erika. camachoOasu. edu 



our motivation we analytically and numerically investigate the mechanism by which solutions 
in certain systems with this resetting conditions stabilize. We further utilize the sociological 
college drinking model as an analogue to analyze certain one-dimensional and two-dimensional 
nonlinear systems, as we attempt to generalize our results to higher dimensions. 

Keywords: College Drinking Model, Logistic Equation, Gompertz Model, Contraction Mapping. 



Introduction 



In many dynamical systems, especially systems found in the biological and natural sciences, it is 
common to find multiple time scales embedded in the same system all of which might be continuous 
or include a combination of both discrete and continuous time. Dynamical systems with multiple 
time scales are also called singularly perturbed systems and slow-fast systems. An excellent example 
of a slow- fast type is the forced Van der Pol equation [5], and is the first system in which the 
phenomenon of chaos was first observed. However for our purposes and interests, we restrict only 
to discrete-continuous time scaled systems. Examples of these arise mostly in ecology. To illustrate, 
recall the famous system of equations known as the Lotka-Volterra model [1,3]. In this model, 
x(t) and y(t) denote, the prey and predator populations, respectively, at any time t. Now let us 
imagine the following prey- predator relationship: The predators are sharks, and the prey arc small 
fish. If there were no more small fish, one might expect that the population of sharks would decline 
in quantity, according to the differential equation 

5 = -cy, c>0. (1) 
at 

When small fish are present in the habitat or environment, however, it seems more reasonable 
that the number of encounters between these two species per unit time is jointly proportional to 
their populations x and y, or in other words, proportional to the product xy, by the law of mass 
action [4]. The addition of this last rate yields an equation for the shark population, and so in a 
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dx 

similar fashion, the dynamics of the prey population — is derived, giving us the nonlinear system 

dt 

of differential equations known as the prey-predator model (2), where a, b, c and d are all positive 
constants. 

dx , 
— = ax — bxy 

'^^ (2) 
dy , 

Although the explicit solution of system (2) is not known, it can be analyzed quantitatively and 
qualitatively [3], that is we can establish some information about the behavior of its solutions, by 
seeking a relation for x and y ignoring temporal space ( a better explanation of this can be found 
in [1] ) . Moreover, we can look at this model in a slightly different way by making the growth rate 
of the prey time dependent, in other words, not constant. If we change the prey population slightly 
as follows, 

^ = a{t)x-bxy, (3) 

then the population dynamics still evolve continuously in time, but what differs now is the birthrate 
of the prey, denoted by some given a{t) which evolves discretely in time. This is due to normal 
ecological changes in temperature or the environment, and as a result the population behaves 
periodically according to this birthrate. In addition, we also find these type of phenomena in fishing 
models with harvesting such as the one found in [1] given by the generalized equation 

f = m-H(t). (4) 

Here H{t) is a non-constant rate of members ( fish ) harvested per unit time, and as in the Lotka- 
Volterra model, the population of fish evolves continuously, yet due to the harvesting rate H{t), 
the population resets periodically, and this periodic reseting is dependent, more realistically, on 
harvesting seasons. 

As observed, such examples of multi-scaled systems are naturally found in nature, but such sys- 
tems are also found in sociological models like a deterministic Model of College Drinking Patterns 
analyzed by Camacho, et al. (2006) [6] ( See also [10] and references therein ). This model has 
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become the onset point of this prospectus journey. A brief description of the continuous determin- 
istic model given by equations (5)-(8) below, is as follows. First a four equation system is given, 
inclusive of an embedded discrete difference equation Rj(tk+i) (9) which is originally called the 
resetting equation due to Scribner, et al. [10]. This equation is crucial to describe periodicity and 
resetting of a population as it takes care of replacing the total number of students that left at the 
end of each academic year, back into the system at the beginning of each academic year by keeping 
the total population in each class constant. Through this equation, every year current variable 
values are reset, resulting in a non-smooth set of solution curves. In other words, the population 
in each class ( TVi, A^2, -^Sj ^4 ) gets reset to bring the total population down to the initial data 
A''(0). This resetting force leads to stabilization but at a different level than the equilibria from the 
continuous system absent from resetting forces. 
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(10) 
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Prom the resetting equation, tk+i and represent the beginning of year fc + 1 and the end of 
year k, respectively, and the negative superscript indicates that tk+i is dependent on the previous 
initial time tk- It can also be noted that N(0) is the total initial population, and cj = ^'-'^ 
j = {1,2,3,4} is a percentage of the total population which enters each class. Moreover, in the 
modified model discussed by Camacho et al. [6], numerical tools were the only option available to 
analyze the nonlinearity of the system, hence numerical simulations have shown that stabilization 
of the solutions occurred even when the continuous system without resetting did not stabilize ( i.e. 
solutions of the system with resetting present did not attain the stable equilibrium belonging to 
the continuous system absent of resetting ). It was also noticed that for certain initial conditions 
the system ( equations (5)-(8) ) stabilized at some artificial equilibrium. Analj^tically, it was 
almost impossible to investigate the behavior in the system as a whole, and the worse was finding 
an explicit solution. This facts have motivated us to search for analytic mechanisms which lead 
solutions of nonlinear systems with different time scales to stabilize periodically in time. In this 
quest, many questions have arised. (?) Under what conditions will the system stabilize and {ii) 
is there a way of proving analytically that stabilization of the system occurs? In the latter case, 
due to the nonlinearities of the four dimensional alcohol model, proving and obtaining conditions 
for stabilization analj^ically are not easy tasks. Thus in this project we have focused on searching 
for mechanisms which stabilize system dynamics, from the simplest model available moving on to 
more complicated systems of larger dimensions once some conclusions about simpler models are 
established. Thus, we begin investigating this mechanisms in nonlinear models case by case. 



Analytical Investigation of Simplest Nonlinear Models 
(Summary of Results) 

In order to describe a mechanism of stabilization in nonlinear dynamical systems, we need to 
start from the simplest linear and nonlinear models available. In this prospectus summary we 
have considered elementary nonlinear cases both one-dimensional and two-dimensional. With the 
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following introductory but significant examples, we show that specific resetting conditions and 
parameter values give rise to the stabilization of these systems. The simplest one-dimensional 

available case considered are exponential growth and decay (Malthus Law) where 

dx , dx 

= ax and — = —ax 
dt dt 

respectively, are the well known differential equations describing exponential growth and decay, 
with their very well known solutions namely 

x{t) = a;(0)e"* and x{t) = a;(0)e-"*. 

Secondly we have considered both, a logistic model 

and a Gompertz (1825) model (taken from [1]) 

^ = -axln{x). (12) 

For these models we have artificially created a resetting condition given by 

a;(T„+i) = 7x{T-+,) (13) 

where 7 is a resetting parameter, and < 7 < 1. On the other hand, in the two-dimensional case 
we have considered a Logistic Coupled system [3] 

= ax(l-^) (14) 



dt \ P 

ft = (1^) 
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and a Gompertz-type System [1] 



dx 



= —axhi {x) 



(16) 



dy 



= — /Syln [x) 



(17) 



dt 



for which we have chosen the resetting condition 



x{Tn+l) 



I 



A 




2/(^n+l) 



(18) 



y(T'n+i) 
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(A more detailed explanation of the mechanisms chosen is given in the proceeding parts of this 
prospectus summary). 

Summciry 

Our illustrations have shown that the resetting mechanisms we investigate, certainly stabilize non- 
linear system dynamics. As a hypothesis we propose the generalization, that any nonlinear system 
can be analytically investigated by following simply the theory of dynamical systems along with 
an advanced calculus approach as follows. It has been established that for a nonlinear autonomous 

system of ordinary differential equations 



X = f(x) 



which has a given initial condition. 



x(0) = xo 
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the fundamental Existence-Uniqueness Theorem [4] states that as long as there exists an open 
subset S of R", which contains xq along with the requirement f e C^{S), then there must be 
some a > such that the initial value problem (IVP) stated has locally an explicit analytical, 
unique solution x(t) on some interval [—a, a]. In addition to the uniqueness of the solution of 
the system, resetting implies intersection of periodic orbits in the state space of the continuous 
dynamical system, with certain lower dimensional subspace called Poincare Sections, transverse to 
the flow of the system. Initial data are thus periodic orbits with initial conditions on the Poincare 
section thereby formulating a Poincare map 

Xn+l = f(x„). 

This iteration process or map allows us to further, invoke the contraction mapping theorem, and 
thus we can claim that the dynamics of the system are guaranteed to stabilize, when a generated 
sequence ( Poincare Sections ) of initial data converges to a fixed point of the corresponding discrete 
mapping. But in general, for n-dimensional systems, it is very diSicult to find a domain of initial 
data, where the discrete map might converge analytically. Investigation of the "size" of this region 
is a typical inverse problem which will be investigated numerically in the following sections. 
Furthermore, there are many possibilities and applications to this work. Higher dimensional sys- 
tems for example, could be analyzed with monotone systems theory, by means of explaining the 
behavior of high dimensional systems by that of one dimensional flow. But in general we want to 
flnd how the interplay of initial conditions and resetting deflnitions affect the behavior of solutions 
in the long run. Ultimately we want to flnd the set of initial conditions that leads to stabilization 
of dynamics by means of a convergent discrete mapping or Poincare Mapping. Our future consider- 
ations include the investigation of any type of model in which resetting occurs and ideally we want 
to find applications to ecological and sociological models. It is also of potential interest to consider 
models in which the periodic administration of a drug, such as chemotherapy for instance, is given 
to cancer patients. 
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Appendix: Examples of Numerical Results 
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Figure 1: College Drinking Model. In this figure, the smooth curves represent solutions of the 
continuous system ( equations (5)-(8) ) without resetting and the the non-smooth curve represent 
the solution of the system with resetting. Note that the stabilization of the continuous system with 
resetting occurs at a level different from the equilibria of the continuous system without resetting. 
That is, when resetting is introduced the system tends to stabilize to an equilibria of its own. 
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Figure 2: One— Dimensional Logistic Model with Resetting Using the resetting condition 
chosen for the one-dimensional logistic model ( equation (11) ), numerical simulations show the 
stabihzation expected, which match those observed by Caniacho et al. Here again, the smooth 
curves correspond to the continuous system without resetting. At some time, the resetting rule 
comes into effect, and the solutions become nonsmooth. Notice in the first case for a resetting 
parameter 7 = 0.67 < 1 and initial condition xq = 0.5 { < xq < 0.9, where a;*^^^;^ — 0.9, and 
^unst — 0), the dynamics of the system with resetting stabilize between the equilibrium of the 
continuous system, yet they follow an equilibrium level of its own. In the second case for initial 
conditions above stable equilibrium, we see the same phenomena ocurring. It is very interesting to 
note that through this resetting mechanism ( or rule ) we have chosen for the simple one-dimensional 
model, stabilization is not the level of the system without resetting. 
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Figure 3: Phase Portrait of Logistic Coupled system absent from resetting The chosen coupled 
system i, = ax(\ — x/P), y — f3xy for < a, /3 < 1 which has explicit continuous solutions given 

by x{t) = —r tR^^^\n.^^ ^^'^ ^ y(0)(l + 2;(0)(e"' - 1)'^^") has a phase portrait where 

equilibria {x* ,y*) = (1, 0), and (0, y*) are a saddle point and a set of degenerate isolated unstable equilibria, 
respectively. Every trajectory is unstable, as the behavior of the system without resetting exhibits a saddle 
node, unless the initial conditions fall on the stable manifold (1,0). 
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Logistic 20 resetting 2{t )=A" z(t" ) with A= [0.5.0.0.5.0,005] 



-X w/ Resetting 

- y w/ resseting 

- X No resetting 

- y No rGEGtting 



Logistic 2D resetting 2(t )=A' z{r ) with A= [1,5,0,0,05.0,005] 
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Figure 4: Two— Dimensional Logistic Coupled System Simulations belonging to the two- 
dimensional system ( equations 14-15 ) again show, with the chosen resetting rule once again the 
system stabilizes, as expected for certain initial conditions and resetting parameter 7. The solution 
y{C) in the continuous system without resetting ( represented by the pink dashed curve ), is tending 
to infinity, while with resetting ( represented by the red solid curve) ?/(t) tends to some nonzero 
equilibrium ( very close to zero ) but in an oscillatory manner with very small amplitude. The 
solution ( represented by the blue solid curve ) in the system with resetting also stabilizes 
to a different level than the continuous solution x(t^ in the system without resetting. Notice the 
amplitude of the oscillations of the solutions of surpass those of y(t) ( the red curve ) and this 
is due to the choice in size of the parameters. It is quite surprising that although the continuous 
system has no stable solution the introduction of a periodic resetting still causes the system to 
stabilize. 
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